InterrillDetachmentRate Subroutine

public subroutine InterrillDetachmentRate(pRate)

compute sediment detachment rate (kg/s). Interril area detachment is considered to be due to the raindrop impact, which breaks the cohesive bonds between particles, making them available for transport.

References:

Foster, G. R., Modeling the erosion process, in Hydrologic Modeling of Small Watersheds, edited by C. T. Haan, H. P. Johnson, and D. L. Brakensiek, pp. 297–383, Am. Soc. of Agric. Eng., St. Joseph, Minn., 1982

Arguments

Type IntentOptional Attributes Name
type(grid_real), intent(in) :: pRate

precipitation rate (m/s)

Variables

Type Visibility Attributes Name Initial
integer, public :: i
integer, public :: j

Source Code

SUBROUTINE InterrillDetachmentRate &
!
(pRate)

IMPLICIT NONE

!Argument with intent(in):
TYPE (grid_real), INTENT(IN) :: pRate !! precipitation rate (m/s)

!Local declarations:
INTEGER :: i,j

!------------end of declaration------------------------------------------------

DO j = 1, interrillErosion % jdim
  DO i = 1, interrillErosion % idim
    IF (interrillErosion % mat (i,j) /= interrillErosion % nodata) THEN             
                                      !conversion m/s -> mm/h 
      interrillErosion % mat (i,j) = 0.0138 * (pRate % mat(i,j) * 3600000. )**2.0 * &
                            rusleK % mat(i,j) * slopeFactor % mat(i,j) * &
                            rusleC % mat(i,j) / 3600. * CellArea (interrillErosion, i,j)
                            
    END IF
  END DO
END DO

END SUBROUTINE InterrillDetachmentRate